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We summarise recent results for the chiral Random Two-Matrix Theory 
constructed to describe QCD in the epsilon-regime with imaginary chemical 
potential. The virtue of this theory is that unquenched Lattice simulations 
can be used to determine both low energy constants S and F in the leading 
order chiral Lagrangian, due to their respective coupling to quark mass and 
chemical potential. We briefly recall the analytic formulas for all density 
and individual eigenvalue correlations and then illustrate them in detail 
in the simplest, quenched case with imaginary isospin chemical potential. 
Some peculiarities are pointed out for this example: i) the factorisation of 
density and individual eigenvalue correlation functions for large chemical 
potential and ii) the factorisation of the non-Gaussian weight function of 
bi-orthogonal polynomials into Gaussian weights with ordinary orthogonal 
polynomials. 

PACS numbers: 02.10.Yn, 12.38.Gc 

1. Introduction 

Non-Hermitian Random Matrix Theory (RMT) has received a lot of 
interest in the past few years due to its relation to QCD with chemical 
potential see [1] a for recent review (and [2j this workshop). Many re- 
sults have been obtained, including correlation functions [3llll[5], individual 
eigenvalues [6] or the phase of the Dirac operator [7], and have been suc- 
cessfully compared to QCD lattice data [8l[6l[9]. We have now understood 
that RMT with (or without) chemical potential is equivalent to QCD [TO] in 
the limit of the epsilon-regime of chiral Perturbation Theory (echPT) [11] . 

The virtue of having /i 7^ is that it couples to F to leading order in 
echPT [12]. The downside of 7^ is of course the sign problem, making 
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unquenched simulations very hard. It was therefore proposed in [13j to use 
imaginary fi instead to determine F, keeping the Dirac operator eigenval- 
ues real and thus making unquenched simulations possible. First results for 
the two-point density [13] were derived directly from echPT for imaginary 
isopin chemical potential and compared to quenched and unquenched Lat- 
tice data (of course real isopin chemical potential could also be simulated 
unquenched). This inspired us to write down and solve the corresponding 
two-Matrix Theory (2RMT) [T3], where in addition partial quenching is 
possible by setting one of the two to zero. This method has already been 
successfully compared to the lattice QCD in [131 115]. In [HJi all unquenched 
density correlation functions were computed (including those for the non- 
chiral theory for QCD in three dimensions). This lead to the construction 
of individual eigenvalues as well [16], and has subsequently been proven to 
be equivalent to the corresponding echPT for all correlation functions |lUj . 

The purpose of this paper is to illustrate the mathematical structure 
of these results by using the simplest possible setting, the quenched theory 
with imaginary /i of isopin type. For results in full generality, including 
partially quenched and unquenched examples we refer to \14\ [T6] . 

The quenched case furthermore helps to point out the differences be- 
tween correlations of two anti-Hermitian Dirac operators with imaginary /i 
isospin and eigenvalues on M?, and one non-Hermitian Dirac operator with 
real /i and eigenvalues on C. Below we show that in the limit of large imag- 
inary fi all correlation functions factorise and become ii-independent, given 
by the product of two single, uncoupled Dirac operators. For large real /i 
however, the complex eigenvalue densities stays ^-dependent and becomes 
rotationally invariant around the origin in C [5]. 

This article is organised as follows. In the next section [2] we recall the 
2RMT and its equivalent echPT, as well as the general results for all corre- 
lation functions. In section [3] we then specify these results to the simplest, 
quenched example with imaginary isospin, including two interesting prop- 
erties. First, the factorisation of all correlation functions for large /x is 
derived and illustrated with several figures. Second, the factorisation of the 
non-Gaussian weight on into two Gaussian weights on M is shown. 



We begin by writing down the partition function of echPT given by |12j 



Here F and T, are the Pion decay constant and chiral condensate, respec- 
tively. In the epsilon regime they have as source terms chemical potential 
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through the charge matrix B =diag(//ilAr;^ , /i2lAf2)) and the diagonal mass 
matrix M =diag({mji}, {771^2}), respectively. 

For this theory all eigenvalue correlation functions are known [14^ [T6] 
and are equivalent [TO] to the chiral 2RMT eq. ([2]). Deriving correlation 
functions from echPT one has to add auxiliary fermion-boson pairs to gener- 
ate the corresponding resolvents, and we refer to [TO] for details. In fact the 
RMT-echPT equivalence holds for any number of chemical potentials, but 
only for two different chemical potentials this theory has been solved. From 
now on we set = —fi2 for simplicity and follow the 2RMT framework as 
it is much simpler. The corresponding partition function is defined as 



Zrmt = / d<i>d^ g-iVTr($t$+v,t^) -Q det[P+ + myi] J] det[P_ +m/2]. 

/1=1 /2=2 

(2) 

The two anti-hermitian Dirac matrices P± are given in terms of two com- 
plex, rectangular random matrices $ and ^ of size x (A^ -|- v) 

^(Dt ± ^^vl/t ) ■ 

When rotating to the eigenvalues xj and yj of 'D± the two random matrices 
get coupled, leading to a non-trivial dependence on the unitary rotations. 
Integrating them out we obtain the following non-Gaussian eigenvalue model 
|14j . up to an overall constant, 

N I ,.00 N2 \ 

Zrmt = n / dxidy,{xiy^Y+^ [| (xf + m)^) JJ (y^ + ni%) (4) 

i=i /i=i /2=i / 

X A^({x2})A^({y2})det 



For later convenience we abbreviate the integrand or joint probability dis- 
tribution function by {?/}). 

2.1. Definitions and Results 
If we define the weight function 

A^i N2 

w{x,y) = {xyy+^ Ylix^ + m}^) Hiy^ + mj^) 

fl=l /2=1 
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we can find the corresponding bi-orthogonal polynomials 

POO 

/ dxdy w{x,y) P„(x^) Qfc(y^) = KSnk , (6) 

^0 

that depend parametrically on the masses. All correlation functions defined 
in eqs. ([8]), ([9]) and ()12p below can then be expressed in terms of the 4 kernels 
K]\f, Hp^, Hn and that are constructed respectively from the two bi- 
orthogonal polynomials Pk and Qk, the polynomials Pk and the generalised 
Bessel transform of its partner 

/•oo 

Xk{x)= dy w{x,y)Qk{y'^) , (7) 
Jo 

the polynomial Qk and its partners transform, and both transforms. The 
density correlation functions are then given by [H] 



(8) 



det 



i=k+l j=l+l 

Hisi{xii,Xi^) MNixi^^yj^) -w{xi^,yj2 



The simplest nontrivial example is the density Ri i{x,y) to find an eigen- 
value of P+ at X and of P_ at y. When all eigenvalues of one kind are inte- 
grated out one finds back the densities of the one-Matrix Theory (IRMT), 
which are then /i-independent. 

Alternatively to the density correlations one can define the so-called gap 
probability that the interval [0, s] is occupied by k eigenvalues and [s, oo) 
by {N — k) eigenvalues of T>^, and that the interval [0,t] is occupied by / 
eigenvalues and [t, oo) by (N — I) eigenvalues of P_: 

rs POO 

Ek,i{s,t) = — — — — — — / dxi...dxk / dxk+i---dxN (9) 



Z{N-k)\{N-ir) Jo 

ft f'OO 

X dxi...dxi dyi+i...dyNV{{x},{y}) 
Jo Jt 

N-kN-l rs rt 

= ., ., / dxi...dxk+i i dyi. . .dyi+jRk+i,i+j .{IQ) 

i=o j=o -^0 -^0 

Here we have also given its expansion in terms of density correlations |16| . 
Obviously if all densities are known all gap probabilities follow, and vice 
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versa. Taking the following derivatives of the gap probabilities 

^ ^dsd' ^ ^' -Pk+iA^^t) -Pk+iA^,t) +pk+i,i+i{s,t)) , 

(11) 

then leads to individual eigenvalue distributions defined as 

Pk,i{s,t) = 2"(/c)w) J dxi...dxk-i J dxk+i...dxN (12) 
X / dyi...dyi^i / dyi+i . . . dy^Vixk = s,yi = t) . 



Jo Jt 

Following eq. (jlOp they can be expanded in terms of densities as well, and 
we will use this expansion below (see eq. p^ ). 

3. The Quenched Theory: Illustrations and Peculiarities 

In the following we illustrate the above results with the simplest quenched 
example, taking Ni = N2 = 0. Prom symmetry the bi-orthogonal poly- 
nomials become equal, Pk = Qk, and they are simply given by Laguerre 
polynomials. Their Bessel transforms become the wave functions of the La- 
guerre polynomials, where details are given in the next subsection 13.21 In 
particular two of the kernels then coincide, Hn = H^. 

The large- iV limit can easily be taken using the standard Bessel asymp- 
totic of Laguerre polynomials. We keep 

= lim 2iV^2 ( = VP'^fij) 

^fc = lim Nxk ( = VT,Xk) , Ck = lim Nyk ( = VT.yk) , (13) 

fixed where in parenthesis the corresponding echPT quantities are given. 
Masses are rescaled as the eigenvalues when present. The limit eq. (jl3p 
results into the following building blocks for the correlation functions, the 
microscopic kernels: 



Jo 



* JumMCt) , (14) 



where k = +1,0,-1 for the limit of Ki\i, and M^v respectively. The 
simplest non-trivial example is the rescaled density Ri^i{x,y): 



/,i,i(e, = x°(e, e)x°(c, o - ea+(e, o x-(c, o 



(15) 
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Pl,l(.S,t) 




pi,i(s = 2,t) eq. (16) 




Fig. 1. The quenched density pi.i(s, t) (top left) vs the individual eigenvalue distri- 
bution pi,i(s, t) (top right) expanded to the order given in eq. pB]) . both ai ly — 
and a = 0.159. The fact that Pi,i(s,t) becomes negative is an artefact of the ap- 
proximation. Higher order terms will ensure that it remains zero at large distance 
from the origin. The lower plots show corresponding 2D cuts at fixed s = 2. The 
advantage oi pi^i(s,t) is that it is localised without background. 



For the corresponding individual eigenvalue distribution we use the expan- 
sion following from eq. (jlOp 



Pi,i{s,t) = />i,i(s,t) - / dx p2,i{x,s,t) - / dy pi^2{s,t,y) + . . . . (16) 

Jo Jo 

Both eqs. p5p and (jlOp are displayed in fig. [TJ For comparison we display 
the same quantities of the /i-independent IRMT. Its rescaled density reads 

piiO = f [J'AO - J.+iiOJ^-iiO] ( = = Pi,o(0 = PoAO ) , 

(17) 

see fig. [2j There we include both the exact distribution of first eigenvalue 

m 



= e-^«' (z^ = 0) , 



(18) 



A07krakow printed on February 2, 2008 



7 




Fig. 2. IRMT a,t v = 0: the quenched density eq. (fTT]) (green), the exact distribu- 
tion of the first eigenvalue eq. (fT8|l (blue), and the expansion [18] corresponding 
to eq. (HH) (red). 



and its corresponding approximation [18] . 



3.1. Factorisation of correlation functions 

In the limit of large chemical potential, a ^ 1, the quenched density 
correlation functions factorise. 



lim Pn,k{Cl,...,n,Cl,...,k) = Y\ CiPn,o{Cl,...,n) TT Ci/50,fc(Cl,...,fc) , (19) 

i 3 

where the two factors are given by the //-independent IRMT quantities 

Pn,0(ei,...,n) = P0.n(6,...,n) = ^j)] . (20) 

This follows from the vanishing of the upper right corner in the determinant 
eq. ([8]) when a is large: the corresponding microscopic kernel M^r converges 
to the weight in this limit 



lim I T-(e,C) - ^lAZ^e-^ 1=0, (21) 



keeping \foL finite. Thus the density from our example eq. (jlSp factorises, 
as is shown in fig. [3] left. The deeper reason for the limit eq. (j2ip will 
become clearer in the next subsection 13.21 

The factorisation of the densities leads to factorised gap probabilities 
and individual eigenvalue distributions as well, as follows from eq. (jlOp : 
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Pi,i{s,t) 




Fig. 3. Factorisation of the quenched density pi.i(s,i): eq. at i/ = and 

a = 4.318 (top left) vs. the factorised eq. ([HI), using eq. fTTl) (bottom left). The 
corresponding approximate individual eigenvalue distribution pi i(s, f) eq. (fT6|) for 
the same values (top right) vs. the factorised exact result eq. using eq. pB]) 
(bottom right). 



E 



(-) 



Ek{s)Ei{t) 



^ I dyi... dyi+j p^j^- (yi , . . . , yj+j) 



(22) 



Differentiating twice as in eq. ()lip we get 



9s df 



k\{pk{s)-pk+i{s)) U{pi{t)-pi+i{t)), (23) 
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and thus from comparing to eq. (jlip 

PkAs,t) = Pkis)pi{t) . (24) 

The IRMT quantities on the right hand side are now exphcitly known, 
without approximations. The comparison in 3D is given in fig. [3] right. 
In order to check the convergence of eq. (fT6]) we can cut the 3D plot and 
compare to the exact factorised result, as shown in fig. HI 




Fig. 4. 2D cut of fig. Otop right a.t s — f .5 (red) vs. the exact factorised distribution 
eq. ^ (blue). 

3.2. Factorisation of the weight function 

In this section we show that the bi-orthogonal polynomials for the non- 
Gaussian weight eq. ([5]) can be constructed in terms of orthogonal polyno- 
mials with Gaussian weight. Our discussion follows closely appendix B of 
|14j . Suppose we have two sets of ordinary orthogonal polynomials 

j dx wi{x) Pkix'^) Pi{x^) = fkSk,i 

j dy W2{y) Qk{y^) Qi{y^) = 9k5k,i , (25) 

with weights wi^2 and norms fk and Qk respectively. Then it follows that 
these polynomials are bi-orthogonal with respect to the weight 

CO ^ 

w{x,y) = wi{x) W2{y) ^ Pk{x ^) Qk{y^) . (26) 

k=o ^^^^ 

In our quenched case we can simply choose the ordinary Laguerre weight 
and its polynomials, 

N 2 

wi{x) = W2ix) = x'^^^^e^'^'^ 

Pk{x') = Qkix') ~ (if;?"') • (27) 
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The identity that allows to link the Laguerre weight and polynomials to the 
non-Gaussian weight eq. ^ is given by 

JVfl+u^) , 2 , 2\ / 1 /i2 

2^2 



w{x,y) = {xyr+' e ^^'^ ' I, [-^Nxy ) (28) 
^0 (n + z.)! " Vl + ^' y VI + ^ 



satisfying eq. (|26p . The Bessel transforms eq. ([7]) then simply result into 
the wave functions 

Xk{y) = j dx w{x,y) Pk{x^) ~ y2.+ig-A^T:^?/^^H J^_^y2^ (29) 

All the kernels can now be easily written in terms of Laguerre polynomials 
and their norms, and we refer to [13] for details. Finally let us reconsider 
the upper left block in eq. ([8]). Due to the above identity eq. ([28]) we obtain 



w{x,y) - MN{x,y) = wi{x) W2{y) Pk{x^) Qkiv^)- (30) 

k=N ^^^^ 

Thus naively taking the limit — > oo we would expect the right hand side 
to vanish. However, in the limit eq. (|13|) this is not the case, and we instead 
obtain M]\f{x,y) — > T~(^,C). Only in the limit of large a ^ 1 limit the 
integral in 2~ extends to oo in the new variables / y/a, making it converge 
to the weight (see eq. ()2ip ). This explains the factorisation in this limit, 
illustrating the subtlety of the large- A^ limit (that is distinguished into weak 
and strong non-Hermiticity for real fi). 
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